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We discuss the structure of parametric resonance which occurs in the process of reheating 
after inflation with two interacting scalar fields. It is found that, for the case of a not too 
large coupling constant, a quasi- homogeneous part of the second, initially subdominant 
scalar field may not be neglected due to its stochastic growth during inflation. This 
fact has strong consequences for the reheating stage: dynamics of background fields 
becomes chaotic after inflation - it consists of subsequent chaotic and regular eras, and 
^ \ a tachyonic instability for inhomogeneous perturbations arises in the quantum reheating 

problem. This instability may pose a problem for the standard reheating scenario. In 
order to avoid it, the coupling constant should be either sufficiently small, or very large. 

O 

^! 1 Introduction 

o 

The simplest versions of the inflationary scenario of the early Universe, which produced definite 
predictions about the total density of matter and the power spectrum of adiabatic perturbations 
in the present-day Universe (later verified by observations), are based on the hypothesis that 
our Universe had passed through a period of quasi-exponential expansion (the de Sitter, or 
inflationary stage) about 14 billion years in the past, before the hot radiation-dominated (RD) 
Friedmann-Robertson- Walker (FRW) era occurred. This inflationary stage was driven by (at 
least) one effective scalar field (named inflaton) which was evolving sufficiently slow during 
inflation. Quantum vacuum fluctuations of this field induced by space-time curvature during 
the de Sitter stage are the source of spatial inhomogeneities of space-time metric and matter 
density which have led to formation of all structures in the present Universe including compact 
objects. 

However, to provide required smallness of initial perturbations, an inflaton field should be 
very weakly coupled to other quantum fields including "usual" ones (i.e., those of the standard 
model of strong, electromagnetic and weak interactions). Also, immediately after the end of an 
inflationary stage, all energy density in the Universe is concentrated in the quasi-homogeneous 
"cold" inflaton field 0. So, to make matching with the subsequent hot RD stage possible, it is 
necessary, first, to transmit the energy density to other quantum fields of matter and, second, to 
heat matter in the Universe. To achieve the first aim, the inflaton should be somehow coupled 
to other quantum fields. The simplest way to model a coupling to bosonic fields of "usual" 
matter is to introduce a second scalar field £ interacting with the inflaton 0. Of course, fermion 
fields and boson fields of a higher spin have to be added, too. Further, we assume the quartic 
interaction Vi nt = ^-</> 2 £ 2 . 
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In the perturbative regime, when occupation numbers of inhomogeneous modes of all fields 
are small, energy transfer from the inflaton to other fields occurs through the process (f)(f) — > 
which is rather slow. The same process, as well the process —>■ (fxp (in the case of the 
-j- inflaton self-interaction) and other non-linear processes, leads to reheating of matter, too. 
However, it was found in 90th that sometimes the energy transfer and reheating may proceed 
much quicker if they occur in a regime of the broad parametric resonance, see u~% (a narrow 
parametric resonance is usually not effective due to expansion of the Universe). The structure 
of the resonance is rather complicated i*: there is a stability - instability chart in the variables 
f 3 ^, kj where k is the inverse wave-length of a £ - mode. 

In the present work, we will show that this picture actually takes place for k = ^- 3> 1 
only. As we will see, in the opposite case the conjecture of a single field driven inflation is not 
correct. This results in a drastic change of dynamics during a reheating stage. Dynamical chaos 
arises in the quasi-homogeneous background problem, that is a general feature for Hamiltonian 
systems with a large number of degrees of freedom (see, e.g., W& regarding chaotic dynamics 
in isotropic cosmology). There is even more dramatic change in the quantum inhomogeneous 
reheating problem: a tachyonic instability arises for sufficiently large K and long wavelengths. 
It may present a problem for the standard reheating model because of generation of strong 
inhomogeneities. To exclude a tachyonic instability without qualitative change of the theory, 
one need the coupling constant k to be small as compared with unity. 



2 Inflationary stage 

Our model of the isotropic flat Universe is similar to that introduced in &: 

4> + + Xcf) 3 + g 2 i 2 4> = , (1) 

£ + 3^ + ^ = , (2) 

#=^(&+ 1 7 i e+fr+&e) • (3) 



The usual point of view EP is that the inflationary stage is driven by the field <p and the 
dynamics of field £ is not important until the inflation is finished. Actually a criterion of 
validity of this statement is given by a value of the effective mass of the field £. If initial 
conditions for (|l] - 0) are such that the mass of £ is much larger than the Hubble parameter 
H, inflation becomes single-field driven soon after its beginning, and the scenario described 
in & is realized. In fact, this condition limits the coupling constant k = g 2 /\: if the inequality 
m \ = 9 2< P 2 ^ H 2 = 27rGA0 4 /3 remains valid until the end of inflation (when Gcf) 2 ~ 1), then 
k> 1 (M 2 =G- 1 ). 

What happens when the coupling constant k is not large? Let us suppose that £ is never 
dynamically important in this regime, too: 

g 2 e « A0 2 . (4) 



Here <p an d £ may be understood both as operators and stochastic quasi-classical quantities. 
Furthermore, let us consider the period of inflation when fluctuations in (f) may be neglected, 
so that the background field <fi may be considered as a deterministic quantity. Its evolution as 
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a function of r = ln(a(t) / a f) , where a/ is the value of the scale factor at the end of inflation 
(thus, t < 0), is given by the equation 

7 M»l. (5) 



TlG 



Then the evolution of £ may be calculated in the Gaussian approximation (i.e., assuming that 
its distribution is Gaussian). The following inequality is valid in the regime under consideration: 
m \ — 9 2 4> 2 ^ H 2 = 27rGA0 4 /3. Since G(f> 2 ^> 1 during inflation, it is satisfied, if g 2 is not too 
much larger than A (in fact, this regime takes place everywhere on the "phase diagram" of our 
theory except for the region k 1!). The equation for the evolution of z = nG((, 2 ) reads: 

dz 2m 2 z GH 2 g 2 z At 2 

— = — I = - 1 . (6) 

dr m 2 4tt At 6tt 2 1 ; 

Let us take z = for |r| = r as an initial condition. Then 

c\t\ 9 ' 2/x — ^P-^, c= Xt ; a n . (?) 

As we can see, there is a stochastic growth of (£ 2 ) due to the second term on the right-hand 
side. Thus, we come to an inconsistency with the conjecture (f|): it is impossible to neglect 
the influence of £ on the large-scale inflationary dynamics, and one needs to solve a two-field 
problem if k is not large as compared with unity. This fact results in a radical change of 
reheating dynamics, too. Now we are going to discuss this question in detail. 
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3 Reheating: solution of a homogeneous problem 

It is convenient & to introduce new "conformal" variables = a<j), £ = a£, r\ = f — soon after 
the end of inflation. The equations ([j] - |3|) take the following form in these variables (a prime 
denotes the derivative with respect to the conformal time rj): 



a" ~ 



__0 + A ^ + ^ = O, (8) 



a 



a" r 



c - -e + 9 2 <rt = o , (9) 

a 



After the end of the inflationary stage, terms containing ^ and ^ may be neglected. Then the 
system (|^-|T0|) is drastically simplified making the problem almost Hamiltonian: 

4>" + \4> 3 + g 2 Z 2 4> = , (11) 

i" + g 2 4> 2 i=0, (12) 

A 2 _ 87T (l (~ A 2 f 1 f ~^2 f \ 2 

I (X 



3Ml \ 2 



+ -(i'T + -4> 4 + g -4> 2 i 2 ) = -^E. (is) 

2^^ T 4 r 2^^ j 3M| v ; 
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Structure of its solutions can easily be understood using the projection of its phase space 
onto the plain 0' = 0, £' = 0. The region where the system moves is bounded by the curve 



(14) 

There are two essentially different regimes of the system behaviour. The first is realized 
when £ and are of the same order. Then both fields oscillate. The system moves somewhere in 
the region near the origin, and its motion is substantially chaotic because of chaotic character 
of energy redistribution between and £ degrees of freedom. For larger k, chaos becomes 
stronger. We call these periods "chaotic eras". 

The second regime is realized when a strong fluctuation of (J^'^J occurs, and the system leaves 
the central region percolating into a "tube" under or below this region. This tube is narrow, 
so the system cannot live inside it during an infinite time. The lower the energy of the system 
E and the higher the control parameter k are, the larger is the number of oscillations inside 
the tube. We call these periods "regular" (because of the possibility of finding an approximate 
solution for and £) or "double - period eras" (because of the analogy with intermittency in 
hydrodynamics) . 
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4 Special solution with synchronous oscillations 

Let us consider a particular case when initial conditions are chosen in such a way that ho- 
mogeneous components of and £ are proportional to each other at the end of inflation. Of 
course, generically it is not so. However, we will see below that this special solution correctly 
reproduces almost all features characteristic for a generic case. 

Examining the equations ( |TT| , |T2"D it is easy to find that if = a£, then the constant a is 
fixed: 

2 9 K / 1C \ 

a = - — r = 7> 15 

g — a k — 1 

and we have the following equation for 0: 

0" + 2 2 03 = o (16) 

Of course, there is no chaos for this solution. Now we consider fluctuations of and £ on this 
background. Introducing the Fourier mode decomposition as usual, we get: 

5i k " + (k 2 + g 2 4> 2 ) 8£ k + 2g 2 ^4> k = , (17) 
Y + (k 2 + 3A0 2 + g 2 i 2 ) 54> k + 2g 2 ^ k = . (18) 



When 5<f) k and 5^ k modes are strongly coupled to each other, it is not clear how to determine 
the number of particles as well as Floquet indices. Nevertheless, it is possible to introduce new 
decoupled perturbations: 

E k = A5i k + B6$ k , (19) 

A fe = C5£ k + D54> k (20) 

where A, B, C, D are some constants to be determined. Omitting a simple calculation, we 
present final equations for these rotated perturbations: 

S^' + (A; 2 + 3^0 2 )S fc = O, (21) 
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K + (k 2 + —g 2 <f> 2 )A k = . (22) 

K 

Therefore, first, the presence of £ 7^ results in rotation in the space of field perturbations. 
Now correct variables are the rotated ones. It should be noted that the definition of fields 
Afc strongly depends on k, so there is no hope to construct something like a stability-instability 
chart for arbitrary initial fields and £. 

Second, the field can be interpreted as an "inflaton" with the wave vector k, because the 
structure of (|2lD is the same as for the inflaton mode if £ = 0. 

Third, there is a tachyonic instability for when n > 2. It means that there is an 
avalanche-like creation of excitations for this mode, but these excitations are quite unlike usual 
particles. Instead, we have a chaotic growth of spatial inhomogeneities. It is interesting to note 
that this special toy solution gives an opportunity to investigate the system behaviour near the 
boundary of chaotic instability. The boundary for the instability for this toy solution is sharp 
- there is no instability for all modes k if K < 2. On the other hand, k > 1 is necessary for this 
solution to exist at all. 

It turns out that almost all these consequences are valid in a generic case (background fields 
and £ are both nonzero and independent) if we consider the WKB regime of an inhomogeneous 
problem. There is only one difference - it is hardly possible to determine an exact boundary 
of instability in this case. This will be discussed elsewhere. 

5 Discussion 

We see that the structure of resonance in the standard reheating theory ([1] - |3]) can be more 
complicated than the picture presented in S*. The stability - instability chart exists for k ^> 1 
only (i.e., in the broad resonance regime), and there is no possibility of its constructing for not 
too large k because of the appearance of chaos in the homogeneous background problem and 
the necessity to rotate modes. The underlying reason is the stochastic growth of (£ 2 ) during 
inflation. As a result, quasi-homogeneous backgrounds for both scalar fields and £ are non- 
zero by the end and after inflation. Thus, chaos occurs in the homogeneous reheating problem, 
which is a usual situation for Hamiltonian system with a sufficiently large dimension of the 
phase space. Note that the Toda-Brumer necessary criterion for the global dynamical chaos (a 
negative determinant of the matrix ||<9 2 [/(0,£)/<90<9£||) is fulfilled for our model. Generically, 
chaotic mixing in the (0, £) plane rapidly provides equipartition of energy density between the 
inflaton (0) and the other matter field (£). Thus, the first aim of reheating is achieved. 

Evolution of conformally transformed fields and £ consists of consequent chaotic (|0| ~ |£|) 
and regular (|0| <C |£|) eras. This phenomenon is similar to intermittency in the theory of 
nonlinear dynamics. There, it is well known that the number of chaotic eras grows while 
increasing an external parameter of a system (1/k in our case). This analogy leads to a 
conjecture that a critical k c exists such that regular eras disappear for k < k c . 

Dynamics of quantum reheating in the case of nonzero £ (t) differs from the picture i-* even 
without dynamical chaos in the problem ((l|-0). The first new feature is a time-dependent 
rotation in the space of modes £0^ and 8^. This effect is not unexpected: while our work was 
in preparation, papers discussing the same effect have appeared Ej©. Its necessity is clear 
from the physical point of view: the Universe should be filled up with particle-like excitations 
after a reheating stage. It is well known that the term "particle" is correctly defined for an 
oscillator with a variable frequency in the WKB-regime only. The correct WKB-solution may 
not be constructed without mode rotation in our case. 
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The second, more important, feature is the appearance of tachyonic instability for sufficiently 
large k and long wavelengths. This effect is different from that recently found in the theory of 
reheating in the course of spontaneous symmetry breaking after the end of the hybrid inflation 
111)0. xhe latter effect occurs in case the oscillation frequency of background fields immediately 
after inflation is much more than the Hubble parameter H at this moment. In our model, this 
is not the case and there is no spontaneous symmetry breaking. So, initially it was hardly 
possible to expect a vacuum instability with respect to the creation of tachyonic excitations. 

The tachyonic instability in our model is more like an instability in chaotic dynamical 
systems with a positive largest Lyapunov exponent []. The question arises immediately: is 
it a desirable or undesirable property? On one hand, this instability is closely connected 
with chaotic mixing leading to a rapid transfer of energy density from the inflaton to other 
fields. On the other hand, it produces strong spatial inhomogeneities in the Universe (and 
inhomogeneously distributed high-energy particles as was in the case of the broad parametric 
resonance type preheating %%). Special investigation is needed to track the further fate of 
these inhomogeneities. At this moment, we want to state only that, to avoid a strong tachyonic 
instability, the coupling constant k should be either sufficiently small, or very large (in the 
latter case, there is no chaos in the homogeneous reheating problem, and the results of & 
remain valid). 

This research was partially supported by the Russian Foundation for Basic Research, grants 
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the Russian Academy of Sciences. 
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Regarding long wavelength perturbations, a close point of view was expressed in the recent paper © which 
appeared when our paper was prepared for publication. 



